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We study the influence of different metallic substrates on the electron induced heating and heat 
dissipation of single Ceo molecules in the junction of a low temperature scanning tunneling micro- 
scope. The electron current passing through the molecule produces a large amount of heat due to 
electron-phonon coupling, eventually leading to thermal decomposition of the fullerene cage. The 
power for decomposition varies with electron energy and reflects the resonance structure partici- 
pating in the transport. The average value for thermal decomposition of Ceo on Cu(110) amounts 
to 21 fiW, while it is much lower on Pb(lll) (2.9 fiW) and on Au(lll) (1.0 fiW). We ascribe this 
difference to the amount of charge transfer into Ceo upon adsorption on the different surfaces, which 
facilitates molecular vibron quenching by electron-hole pair creation. 



INTRODUCTION 

A single molecule junction that is exposed to the flow 
of an electron current will experience an increase of tem- 
perature due to the heat generated by conduction elec- 
trons. Recent theoretical studies predicted that Joule 
heating in molecular junctions can be large enough to 
affect the reliability of molecular devices pQ. The tem- 
perature of a molecular junction is difficult to estimate. 
Recent experimental approaches have found agreement 
in that a single molecule device can easily reach several 
hundreds of degrees under normal operation (flow of mi- 
crowatts when powered with 1 Volt) [2j 02 [H [5] . For most 
molecules, such temperature reaches the limit of thermal 
degradation. Therefore, a deeper knowledge of the mi- 
croscopic mechanisms of heat generation and dissipation 
on the nanoscale is required in order to improve the per- 
formance of a molecular device. 

The temperature at the junction is a consequence of 
an equilibrium between heating and dissipation of heat 
away from the molecule. Heat generation is essentially 
caused by scattering of electrons with molecular vibra- 
tions. Heat dissipation is expected to follow several 
mechanisms involved in the coupling of the hot molecule 
with the substrate acting as an external cold bath. Ac- 
cording to this, the temperature reached by a single 
molecule during electron transport should depend on the 
lead's material. Furthermore, one would expect that 
both the phonon bandwidth and the electronic density 
of states around the Fermi energy should govern a priori 
the mechanisms of heat dissipation (cooling) and, hence, 
the molecular temperature. 

In this paper we analyze the effect of the lead's material 
on the heat dissipated during electron transport through 
a single molecule. We use a low temperature scanning 
tunneling microscope (STM) to inject an electron cur- 
rent through a single fullerene adsorbed onto clean metal 
surfaces of different nature. Rather than measuring the 
temperature reached by the single molecule, we follow the 
approach described in our previous work [5] . We look for 



certain current and sample bias (V s ) values at which the 
molecule under the STM tip undergoes an irreversible 
thermal degradation. The applied power required to 
reach the decomposition limit (Pdec) varies slightly with 
the electron energy (eV s ) in a sequence that resembles 
the energy alignment of the fullerene resonances for ev- 
ery metal substrate, confirming a resonant mechanism of 
electron heating [6]. However, we also find that the mean 
values of Pd ec (Pdec) depend more strongly on the sub- 
strate onto which the fullerenes are adsorbed. The way 
that Pdec scales for different metals cannot be explained 
in terms of intrinsic properties of each material. Instead, 
the observed trend seems to be a result of the mecha- 
nism of interaction between metal and molecule. In par- 
ticular, we find a correlation between the behavior Pd ec 
and the amount of charge transfer into molecular unoc- 
cupied states, as resolved using local spectroscopy mea- 
surements of Ceo molecules on the different substrates. 
Hence, our results suggest that charge transfer processes 
at the metal-molecule contact provide an effective path- 
way to dissipate heat from a hot molecule through the 
creation of e-h pairs. 



ADSORPTION OF C 60 ON CU(llO), PB(lll) 
AND AU(lll) 

The experiments were performed in a custom-made 
ultra-high vacuum STM, which is operated at a tempera- 
ture of 5 K. We chose Cu(110), Pb(lll) and Au(lll) sin- 
gle crystals as substrates because of their different prop- 
erties regarding chemical reactivity, as well as phonon 
band widths and densities of states at the Fermi energy 
(Ep). Cleaning of the metal substrates was performed 
by standard sputtering-annealing cycles under ultra-high 
vacuum, ensuring an atomically clean and flat metal sur- 
face. Indentations of the STM tip into the substrate while 
applying a tip-sample bias were used to clean the tips 
in regular intervals. The tips are therefore believed to 
be composed of the same material as the substrate. A 
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sub-monolayer coverage of Ceo molecules was deposited 
from a Knudsen cell on the metal surfaces at room tem- 
perature. On Cu(llO) the system was further annealed 
to 470 K to ensure that the molecules self-assemble in 
ordered domains and populate a thermally stable ad- 
sorption state. The structure and electronic configu- 
ration of the molecular layers were characterized using 
STM and scanning tunneling spectroscopy (STS) mea- 
surements (Fig. [T]). 
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FIG. 1: (a-c): STM images of C 60 islands on Cu(110) (h = 
1.0 nA, V s = 2.25 V), Pb(lll) (I t = 2.0 nA, V s = 0.5 V), 
and Au(lll) (It = 0.6 nA, V s = 0.4 V). (d-e): Correspond- 
ing tunnel conductivity spectra (dI/dV s ) of Ceo on the dif- 
ferent surfaces. Bars mark the fitted positions of HOMO, 
LUMO and LUMO+1. Two peaks were fitted to the LUMO 
and LUMO+1 levels on Pb(lll). The spectra were measured 
by positioning the STM tip on top of a single molecule and 
ramping V s while keeping the tip-molecule distance constant 
(feedback- loop open). dI/dV s data were obtained by using a 
lock-in amplifier with a rms modulation amplitude V ac . (Cu: 
Rjunct = 1.1 GVt, Vac = 20 mV, Pb: R junct = 0.7 GVt, 
Vac = 5 mV, Au: R junct = 0.3 GQ, V ac = 30 mV .) 

On Cu(110), Ceo forms ordered islands with a pseudo- 
hexagonal structure (Fig. [I] (a)), in which the fullerenes 
adopt a well-defined adsorption configuration with a 
pentagon-hexagon C-C bond pointing upwards [5j 0. 
An STS spectrum on these molecules (Fig. [I] (d)) re- 
solves a clear spectroscopic fingerprint characterized by 
a sharp resonance at ~1.5 eV above the Fermi level 
and associated to the alignment of the LUMO+1 res- 
onance (LUMO: lowest unoccupied molecular orbital). 
The LUMO resonance appears as a broader peak cen- 
tered at ~0.2 eV and is partially occupied. 

Fullerene adsorption on Pb(lll) at room temperature 
results in highly ordered hexagonal islands (Fig. [l](b)). 
Here, scanning tunneling spectroscopy reveals the LUMO 



and LUMO+1 derived resonances as in Cu(110) but more 
pronounced and with a characteristic split structure due 
to the breaking of degeneracy upon adsorption [8] (Fig. [I] 
(e)). The energetic alignment of the LUMO close to Ep, 
with a small tail crossing it, indicates a small amount of 
charge transfer into Ceo- 

Ceo islands on Au(lll) evidence a similar hexagonal 
lattice as on Pb(lll) (Fig. [l](c)), but with large variety 
regarding molecular orientations [9 , as one can deter- 
mine from intramolecular structure resembling the lobed 
shape of the LUMO resonance [10]. Both LUMO and 
LUMO+1 resonances are resolved in STS spectra as pro- 
nounced peaks, independently of the molecular orienta- 
tion [TT] , The LUMO peak appears located far from Ep, 
at ~ 0.8 eV, thus indicating that charge transfer from 
the substrate is very small [T2] . 



C 60 DECOMPOSITION ON METAL SURFACES 

In our experiment, we approach the STM tip a distance 
Z towards a single Ceo molecule while holding the sample 
bias V s constant. During approach we record the current 
flowing through the molecular junction (I(Z)). The tun- 
nel regime is revealed by an exponential increase of I(Z) 
with diminishing gap distance. At a certain approach 
distance the I(Z) curves deviate smoothly from the expo- 
nential dependence, indicating the onset of a tip- molecule 
contact [3] . The conductivity at this point is a small frac- 
tion of Go (Go = 77.5 nS). For small positive sample bias 
the molecule remains intact upon contact formation and 
even further indentations of several Angstroms, leading 
to a stable junction with the molecule contacted by the 
STM tip on the one side, and the metal surface, on the 
other [3[5]. In this case, the integrity of the indented 
molecule after tip contact can be verified by its appear- 
ance in the STM images and, especially, by its electronic 
fingerprint in dI/dV s spectra. 

Fig. [3] shows, for the case of Pb(lll), the effect of ap- 
proaching a molecule at bias voltages above a threshold 
value of V s = 1.0 V. A sharp drop appears in the I(Z) 
curves before reaching the contact regime [13], denoting 
the occurrence of an irreversible change in the junction 
(Fig. [3](c)). Similar results can also be found on Cu(110) 
for V 8 > 1.2 V [5, and on Au(lll) V s > 1.5 V. After 
this current drop STM images reveal that the molecule 
is transformed into a lower feature (Fig. [3] (d)) and the 
characteristic resonances of the fullerene icosahedral cav- 
ity are absent from its corresponding STS spectrum (Fig. 
[3] (e)). Hence, the discontinuous current drop is a finger- 
print of degradation of the fullerene cage. The precise 
way in which the Ceo decomposes can not be determined 
in detail in our measurements. The effect is observed 
solely on the molecule selected for the tip approach and 
shows a high reproducibility, as depicted in Fig. [3] (f) 
for the case of Ceo on Pb(lll) and as shown in Ref. [5], 
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FIG. 2: Conductance and current vs. Z-distance plot I(Z) on 
top of a C 60 molecule on Cu(llO) (V s = 0.5 V). The blue 
shaded area indicates the contact regime. 



on Cu(110) [13] . A mechanical breaking of the fullerene 
cage can be excluded since the process takes places be- 
fore a tip-molecule contact is formed (i.e. in the tunneling 
regime). We can then conclude that the decomposition 
is a current induced thermal process. 

The current Id ec at which the Ceo decomposition oc- 
curs can be taken from the I(Z) approach curves. Id ec 
turns out to show very reproducible values distinctive 
for every metal substrate and a clear dependence on the 
bias voltage used during the tip approach. Fig. [4] shows 
the statistical average of Idec plotted vs. V s for all three 
surfaces. Adsorbed on Cu(110), Ceo molecules can with- 
stand currents of several tens of microamperes, typically 
one order of magnitude larger than for Ceo on Pb(lll) 
or Au(lll). In general, Id ec decreases monotonously as 
the applied bias is increased. Additionally, some faint 
steps can be inferred. This substructure becomes more 
evident when we plot instead the applied power neces- 
sary to decompose the molecule, Pd ec — V s x 7^ ec (Fig. [3] 

(b-d)). 

The monotonous decrease of Idec transforms in a fairly 
flat behaviour of Pd ec with the applied bias. Fig. [3] (b-d) 
shows that the average power applied for degradation de- 
pends strongly on the substrate used: Pdec ~ 21 fiW for 
Cu(110), Pdec ~ 2.9 /iW for Pb(lll), and P dec ~ 1 /iW 
for Au(lll). Superimposed to these values, the stepped 
substructure appears now more pronounced. The struc- 
ture is correlated with the alignment of the unoccupied 
states of the Ceo molecule, also shown in Fig. [3] (b-d). 
In general we find that Pd ec decreases whenever a new 
resonance (here the LUMO+1 and LUMO+2) enters the 
conduction energy window. Hence, the molecular reso- 
nance structure is reflected in the power needed to de- 
compose a molecule. Next, we interprete the origin of 
this resonant substructure (section 4) and of the strong 
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FIG. 3: Moleculare decomposition of Ceo on a Pb(lll) sur- 
face, (a, b): STM image and dI/dV s plot of an intact Ceo 
(central molecule), (c): Typical I(Z) curve of a tip approach 
experiment showing one decomposition event in the tunnelling 
regime (V s = 2.75 V). (d, e): STM image and dI/dV s plot 
of the molecule from (a, b) after decomposition events on the 
central molecule and on molecules at the image border. A 
height difference of 0.7 A and the disappearance of the charac- 
teristic LUMO and LUMO+1 resonances confirms the degra- 
dation of the Ceo molecule. ( (a): It — 0.2 nA, V s = 0.5 V, 
(d): h = 0.2 nA, V s = 2.25 V, (b,e): R jU nc = 2.3 GQ 
) (f): Image of the letters "STM" created by successive 
decomposition events, to illustrate the large reproducibility 
(It 0.2 nA, V s 0.5 V). 

substrate dependence of the applied power for degrada- 
tion Pdec (section 5). 



MOLECULAR HEATING AND COOLING 
MECHANISMS 

The origin of the sub-structure in the Pd ec vs - V s plots 
can be understood from current theoretical models de- 
scribing the electron-induced heating of single molecules 
during electron transport [15j [16]. Heat inside the 
molecule is generated by inelastic scattering of tunnel- 
ing electrons with molecular vibrations (Fig. [5] (a)). For 
the large tunneling rates used in our experiment elec- 
tron scattering leads to a non-equilibrium distribution of 
excited modes, whose internal energy can be associated 
with a molecular temperature T m . For a certain set of 
current and bias values, T m depends on the balance be- 
tween the heat generated by the inelastic scattering of 
electrons with molecular modes, and heat dissipated into 
the "cold" electrode, in our experiment at 5 K [6 , 15] [17]. 
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FIG. 4: (a): Statistical average of decomposition currents 
Idec(Vs) in the tunnelling regime for the three different sur- 
faces. The total numbers of recorded events are: 152 on 
Cu(110), 150 on Pb(lll) and 105 on Au(lll). (b-d): Sta- 
tistical average of the decomposition power Pd ec (V s ) for (b) 
Cu(110), (c) Pb(lll) and (d) Au(lll). 2a error bars are in- 
dicated in all figures. Normalized dl/dV spectra of Ceo on 
the different substrates are added as shaded curves for com- 



parison (Cu: Rjun 
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Recent calculations [6 have shown that when a new 
resonance level E r enters into the transport energy win- 
dow (hence, when eV s > E r ) a steep temperature in- 
crease takes place in the molecule as a consequence of 
more vibronic levels being accessible (Fig. [5] (a)). Hence, 
the increase in temperature upon crossing a resonance 
level manifests itself as a step-like decrease of the corre- 
sponding decomposition current Id ec and power Pdec, as 
shown in the plots of Fig. [4j 

In contrast to this, tunnel electrons can also absorb vi- 
brons of a hot molecule (Fig. [5](b)), leading to a reduced 
heat generation. When the sample bias lies right below 
a resonance level eV s < E r ) this cooling mechanism can 
be very effective and, eventually, causes a lower rise of 
molecular temperature with V s [5] . This behavior trans- 
lates into plateaus in Id ec an d increase in Pd ec a t bias 
values below the corresponding resonance energy. 



(a) 





FIG. 5: Schematical drawings of vibron interactions with tun- 
neling electrons for a molecule inside a tunnel junction: (a) 
Inelastic electron-vibron scattering at a resonance level E r 
heats the molecule; (b) Vibron assisted tunneling is a cooling 
effect most effective for electron energies close to and below 

Qualitatively, the combination of both heating and 



cooling processes accounts well for the correspon- 
dence between the stepped (oscillating) behavior of 
the Idec(Pdec) plots and the resonant structure of the 
molecule in Fig. [4] [18 . However, it cannot explain the 
striking differences of Pdec found for the different sub- 
strates. Since the mechanisms of heating/cooling by tun- 
neling electrons depend on the coupling between tunnel 
electrons and molecular vibrons they are not expected to 
vary much on the different substrates. Hence, the large 
differences in Pd ec must be related to substrate mediated 
mechanisms of cooling the molecule. 



SUBSTRATE DEPENDENCE OF 
DECOMPOSITION POWER 

Further mechanisms of molecular cooling are those 
in which molecular hot modes decay by creating quasi- 
particle excitations in the substrate (either e-h pairs or 
phonons). In agreement with our findings, these excita- 
tions do not depend on the sample bias but may vary de- 
pending on the substrate's electronic and phononic prop- 
erties. 
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FIG. 6: Schematical drawings of mechanisms of molecular 
vibration decay into excitations of the cold metal substrate: 
(a) decay into the metal phonon band and (b) creation of 
excited metal electrons (electron-hole-pairs) close to the Fermi 
edge are substrate dependent cooling effects. 



Decay of molecular vibrons into substrate phonons 
plays a minor role in cooling the molecule, since this 
mechanism is limited by the phonon band width of the 
substrates (Fig. [6] (a)). The Cu(110) phonon band is 
approximately 42 meV wide and therefore larger than 
in Pb(lll) (14 meV [19 ) and in Au(lll)(16 meV [20]). 
However, the 174 internal modes of Ceo have energies be- 
tween 33 meV and 200 meV. Hence, substrate phonons 
are expected to be primarily coupled to external molec- 
ular vibrations of the Ceo molecule with respect to the 
surface, which do not contribute to the thermal decom- 
position of the fullerene cage. 

The most important contribution to the observed sub- 
strate dependence of Pd ec is the varying efficiencies of 
molecular cooling through electron-hole pair creation 
(Fig. [6] (b)) [21 . Since the final states lie at excita- 
tion energies in the order of the vibrational band width 
of the molecule, the decay rate scales with the density 
of states around Ep. On Cu(110) and Au(lll), surface 
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states provide most of the available states in this energy 
region, according to available data of DOS at Ep [22] . 
It can therefore be expected, that e-h pair creation is 
favoured on Cu(llO) with respect to the gold surface, in 
agreement with our observations. However, this mecha- 
nism fails to explain the intermediate Pd ec value found for 
degradation of Ceo on Pb(lll), since this metal surface 
has no surface state, and whose density of bulk states at 
Ep is also the lowest of all three metals [22] . 

To solve this puzzle we note that the absorption of 
molecular vibrons through excitation of electron-hole 
pairs rather depends on the density of states around Ep in 
the Cqo molecule coming from the surface [5j[23]. There- 
fore, this dissipation channel should not be viewed as 
an intrinsic property of the substrate material, but as a 
consequence of the molecule-surface interaction. In our 
case, we could experimentally determine that the LUMO 
derived resonances of Ceo on different metals exhibit a 
different degree of weight at the Fermi energy, associated 
to different amounts of charge transfer from the surface. 
In fact, the degree of charge transfer (increasing from 
Au(lll) over Pb(lll) to Cu(110),_as shown in Fig. [F] 
(a-c)), follows the trend found for P dec in Fig. g 
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FIG. 7: (a-c): Magnification of the dI/dV s curves of '. 
(d-f) showing the LUMO of C 60 on (a) Cu(llO), (b) Pb(l 
and (c) Au(lll). The occupied portions of the LUMO peaks 
are marked as shaded areas. From this a decreasing amount 
of electron transfer into the molecule can be expected for 
Cu(llO) to Pb(lll) to Au(lll). (d): Schematical drawing 
of the proposed mechanism of vibron decay by electron-hole- 
pair creation in a partially occupied molecular resonance. The 
large spacial overlap of molecule and metal states at Ef in- 
creases the decay probability. The e-h-pairs excited inside the 
molecule can leave before they recombine. 

This allows us to depict a qualitative picture, in which 
charge transfer from the metal causes an increase of the 
density of molecular states around E F . This favours the 
generation of e-h pairs at the molecule, which are then 
reflected back into the metal where they recombine. The 
rate of heat dissipation thus increases with the charge 
transfer, causing that a larger power is required to ther- 
mally degrade the single molecule. Ignoring additional 



mechanisms and effects like hybridization of metal and 
molecule states, which probably could also play an im- 
portant role here, our results hint that charging a single 
molecule in contact with a metal electrode can help to 
sustain larger current densities crossing through a single 
molecule. 



SUMMARY 

Our experiments reveal that electronic currents in the 
range of 0.1-20 \i A and powers in the order of 1-30 fiW are 
sufficient to generate heat in single Ceo molecules lead- 
ing to thermal decomposition on metal surfaces. The de- 
composition power results from the balance of molecular 
heating and cooling. While the former is substrate inde- 
pendent, the latter varies on the different substrates as 
a function of charge transfer into the molecule. Charge 
transfer assures an effective quenching of molecular vi- 
brons into electron-hole pairs. In order to increase the 
current density a molecular junction can sustain, it is de- 
sirable to choose a substrate where molecular adsorption 
leads to a partial filling of molecular states by substrate 
electrons. 
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